Prepared for submission to JCAP 



Cosmological perturbations in 
theories with non-minimal coupling 
between curvature and matter 



Orfeu Bertolami, a l Pedro Frazao 6 and Jorge Paramos 6 

a Departamento de Fisica e Astronomia, Faculdade de Ciencias, Universidade do Porto, 
Rua do Campo Alegre 687, 4169-007 Porto, Portugal 

^Instituto de Plasmas e Fusao Nuclear, Instituto Superior Tecnico, Universidade Tecnica de 
Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 

E-mail: orfeu.bertolami@fc.up.pt, pedro.frazao@ist.utl.pt, paramos@ist.edu 

Abstract. In this work, we examine how the presence of a non-minimal coupling between 
spacetime curvature and matter affects the evolution of cosmological perturbations on a ho- 
mogeneous and isotropic Universe, and hence the formation of large-scale structure. This 
framework places constraints on the terms which arise due to the coupling with matter and, 
in particular, on the modified growth of matter density perturbations. We derive approximate 
analytical solutions for the evolution of matter overdensities during the matter dominated era 
and discuss the compatibility of the obtained results with the hypothesis that the late time 
acceleration of the Universe is driven by a non-minimal coupling. 
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1 Introduction 

Over the past decades, the cosmological perturbation theory, resulting from the linearized Ein- 
stein field equations, has become a cornerstone of modern cosmology [1]. It aims to describe 
the large-scale structure of the Universe (LSS) [2] and to match the observed fluctuations 
of the cosmic microwave background (CMB) [3]. In fact, the Universe is, on large scales, 
homogeneous and isotropic since primordial times. However, small primordial density pertur- 
bations, resulting from vacuum fluctuations during inflation, grow as the Universe expands 
via gravitational instability — resulting in the structures one may observe at the present time, 
such as galaxies and clusters of galaxies. These matter density perturbations evolve during 
cosmic expansion along with metric perturbations around this smooth spacetime. 

The late-time accelerated expansion of the Universe [2, 3] can be explained either by an 
exotic form of energy, the so-called dark energy, or by a modification of gravity at cosmological 
scales. In particular, in modified theories of gravity such as f(R) gravity [4], where the linear 
dependence of the Einstein-Hilbert action on the Ricci scalar is substituted by a non-linear 
function, it has been shown that the new terms that arise in the modified field equations can 
effectively lead to late-time cosmic acceleration [5-7]. However, describing the cosmological 
evolution of the accelerated expansion of the background is not sufficient to break down the 
degeneracy between modified theories of gravity and usual dark energy models — an issue 
further complicated by the possibility to link f(R) models with appropriate scalar-tensor 
theories via a conformal transformation [8]. 
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The growth of matter density perturbations depends strongly on the underlying theory 
of gravity (namely General Relativity (GR) or a modification of it), since the modified field 
equations lead to different evolutions of cosmological perturbations and alter the characteristic 
imprints that these leave in the CMB and in the matter power spectrum inferred from galaxy 
clusters. It is therefore an outstanding probe to modified gravity models, allowing for the 
discrimination between the latter and usual dark energy models (which assume the validity 
of GR) [9-14]. For instance, in the context of f{R) theories, it is not possible to describe 
simultaneously small scales of galaxy surveys data and large scales of CMB data, as the 
latter are suppressed in comparison with the former [15]. The cosmic microwave background 
high-precision data can also be used to constrain f(R) gravity models [16-22]. 

By the same token, the evolution of matter perturbations can also lead to stringent 
restrictions on f(R) theories that reproduce the Einstein-Hilbert action at high curvatures 
through a completely general procedure [23-25]: in particular, the linear cosmological evolu- 
tion of f(R) gravity can reduce the large-angle CMB anisotropy, as well as induce a different 
profile for the linear matter power spectrum and the correlations between the Integrated- 
Sachs- Wolfe effect and the angular power spectra of galaxies [23]. 

Another interesting gravity probe is gravitational lensing [26]: indeed, the relation be- 
tween the gravitational potentials in the metric, which are responsible for gravitational lensing, 
and the matter over densities, depends on the theory of gravity [27]. 

In the context of the non-minimal coupling between scalar curvature and matter [28], it 
has been shown in several works that both dark energy and dark matter can be effectively de- 
scribed within this framework, along with other cosmological and astrophysical implications. 
For instance, the mimicking of dark matter [29, 30] or even dark energy [31], the reheating 
scenario after inflation [32], and the change in the gravitational potential [33]. 

In this work, one examines the possibility of a non-minimal coupling between curvature 
and matter via the study of the evolution of cosmological perturbations, by following the 
implications of the modified field equations and the non-conservation of the energy-momentum 
tensor. 

This work is organized as follows: first, one introduces the non-minimal coupling model. 
Then one shortly reviews the formalism of the cosmological perturbation theory and obtains 
the linearized equations in the presence of a non-minimal coupling. In the last sections, 
one derives the equations of motion for scalar perturbations and matter overdensities in the 
subhorizon approximation and studies the respective evolution for the relevant modes in the 
formation of LSS along with the respective conclusions. Finally, conclusions are drawn as to 
the main effects stemming from the presence of a non-minimal coupling. 

2 The non-minimal coupling 
2.1 Modified field equations 

A recent generalization of modified theories of gravity resorts to a non-minimal curvature- 
matter coupling, further extending the presence of non-linear functions of the scalar curvature 
in the action functional [28]. The considered action is 



where fi(R) (with i = 1,2) are arbitrary functions of the scalar curvature and g is the metric 
determinant. Setting f2(R) = 0, one obtains the usual f(R) theories of gravity, and GR is 




(2.1) 
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recovered if in addition the linear function fi(R) = R/k is considered (with k = 8nG), where 
G is the Newtonian gravitational constant. 

The field equations are obtained from the variation of Eq. (2.1) with respect to the 
metric [28], 

FK ~ \sth - (9^V P V U - 5£D) F = (l + f 2 ) T>? , (2.2) 

where T„ is the energy-momentum tensor, = dfi/dR and F = F\ — IFip, with p the 
energy density. One assumes a matter Lagrangian of the form C m = —p (see Refs. [34] for a 
thorough discussion). 

Taking the trace of Eq. (2.2), one gets 



FR-2f 1 + 3DF=(l + f 2 )T , (2.3) 

where T is the trace of the energy-momentum tensor. In the following, one considers that 
the matter-energy content of the Universe is well described by a perfect fluid, T^ u = (p + 
p)u^u u + pg^ u , where is the four-velocity field in comoving coordinates (with u = 1 from 
the condition u^u^ = — 1). 

For an homogeneous and isotropic spacetime, described by a Friedmann- Robertson- 
Walker (FRW) metric, the background evolution of the field equations for the non-minimal 
coupling, Eq. (2.2), takes the form 



3FH 2 = p + p c , (2.4) 
-F(2H + 3H 2 ) =p + Pc , (2.5) 

where the new terms are due to the modification of gravity and appear as curvature compo- 
nents of the energy-density and pressure, which are defined as 

Pc = -~(f-FR)-3HF , (2.6) 
Pc = f2P+(F + 2HF)+ l -{h-FR) , (2.7) 

where the dots refer to time derivatives and f = fi — 2f<iP- 
2.2 Non-conservation of the energy-momentum tensor 

One distinct signature of the non-minimal coupling is the non-conservation of the energy- 
momentum tensor [28]. Using the Bianchi identities, V^G> = 0, and taking the covariant 
derivative of the field equations, Eqs. (2.2), results in 

V^ = -( / 9^ + r^)V M log(l + / 2 ) , (2.8) 

which one may identify as an energy-momentum exchange between matter and spacetime 
geometry. Setting f2(R) = leads to the usual GR covariant conservation of the energy- 
momentum tensor. 
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However, in a perturbed spacetime, these covariant derivatives will differ from that of 
the background metric due to the perturbation on the Christoffel symbols, and one must also 
take into account the perturbations on the r.h.s. of Eq. (2.8), which are non-vanishing due 
to the presence of the non-minimal coupling. 

In this context, for a perfect fluid description, the above non-conservation of the energy- 
momentum tensor, Eq. (2.8), leads to the extra force [28] 



where W = g^ u +u^u u is the projection operator. This is orthogonal to the fluid four- velocity, 
since f^u^ ~ h^ u u^ = 0. 

3 Cosmological perturbations 

The theory of cosmological perturbations resorts to the evolution of small fluctuations around 
an idealized homogeneous and isotropic background Universe. These cosmological perturba- 
tions can be decomposed in three types: scalar, vector and tensorial modes, which have 
different and independent evolutions during cosmic expansion. In general, only scalar per- 
turbations present instabilities that grow during the expansion of the Universe and lead to 
the formation of LSS. Conversely, vector and tensorial modes do not present such growing 
inhomogeneities: the former decay with the cosmological expansion and the latter do not 
present any coupling with the energy density or pressure because the non-diagonal terms of 
the energy-momentum tensor are zero. Hence, one restricts the attention to the scalar mode 
of cosmological perturbations. 

The linearized modified field equations and the non-conservation of the energy- momentum 
tensor constitute the set of equations which govern the evolution of cosmological perturba- 
tions around the smooth FRW background. One now considers the most general form of the 
line element for the perturbations of the FRW spacetime, which in the longitudinal gauge 
yields the line element [1], 



where the perturbations coincide with the Bardeen gauge invariant potentials. 

The linearized field equations, Eq. (2.2), and the non-conservation of the energy- 
momentum tensor, Eq. (2.8) lead, to first order, to a set of equations governing the cos- 
mological dynamics of the potentials $ and along with the matter density and pressure 
perturbations, as one will derive in the following paragraph. For a further discussion on a 
more general coupling function of the matter Lagrangian see Ref. [35]. 

3.1 Linearized field equations 

The linearization of the field equations is obtained considering the first order perturbation 
in the curvature due to metric fluctuations. Therefore, the time-time component of the field 
equations, Eq. (2.2), is 




(2.9) 
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(1 + h) &P - (F 2P ) 5R 



(~2F 2 p) 



-3* - 3H($ + 2*) - 3(H + H 2 ) (2$) - -= V 2 $ + 3{H + H 2 ) 5 In (F 2 p) 



6H(^ + H$) 



_V 2 (2^) + 3(i7 + // 2 )(51iiFi 



+ 



1 



6HF<f> + F($ + 3*) -3H5F + -^V 2 (8F) 



(3.2) 



while the spatial-spatial component reads 



5)\{l + h)5p+p5f 2 } 



-2F 2P ) 



1 



- (H + 3H 1 ) (2$) - H($ + 6*) + -~ 



+ 5!- 



□F + (JET + 3# 2 ) (SF) 



F 

a 



2 5 lk V k [Vj($-^ + (JlnF)] 



(3.3) 



Here, one identifies three types of terms in the r.h.s. of these equations: one due to the 
presence of the coupling function, f 2 (R), other stemming from the fi(R) non-linear curvature 
term, and a third involving a linear combination of these, present in the function, F = 
F\ — 2F 2 p, and respective fluctuation 5F. 

3.2 Matter density perturbations 

After considering the linearization of the modified field equations, one now examines the linear 
approximation of the components of the non-conservation of the energy-momentum tensor, 
Eq. (2.8). First, it is relevant to introduce a new function which may be identified in the 
extra force given in Eq. (2.9) for a pressureless perfect fluid 



/" = - (3„* c ) hP" 
where one defines a coupling potential 



(3.4) 



* c = log(l + / 2 ) 



(3.5) 



i.e. the extra force / M is the gradient of this potential. Scalar fluctuations of the metric will 
lead to a corresponding perturbation in this potential, 5$> c = log' (1 + f 2 ) 5R. 

The equation of motion for matter overdensities results from the linearization of the non- 
conservation of the energy-momentum tensor, Eq. (2.8): the r.h.s. of the time component is 
zero, while the spatial components read 



V M T? = -(5? Pm + T?)d^ c , (3.6) 

where one takes the case of a Universe dominated by matter, p ~ p m . Introducing the velocity 
potential v, defined through U{ = —diV, one may calculate from the above expression the first 
order perturbation in the components 5Tq of the energy-momentum tensor 



- 5 - 



v - $ c v = $ — 8<f> c 



(3.7) 



Thus, the time and spatial components of the non-conservation of the energy-momentum 
tensor, Eq. (2.8) lead to a set of coupled differential equations which determine the evolution 
of the matter fluctuations, 8p m : 

5 P^ + M 5 ^ = M-%v m , (3.8) 

Pm Pm Q> 

v m + (H + § c )v m = + 5$ c ) , (3.9) 
a 

where the last equation is obtained by introducing the velocity potential, v m = v/a, in 
Eq. (3.7). These equations, along with the perturbed field equations, Eqs. (3.2) and (3.3), 
constitute the set of equations that govern the cosmological evolution of scalar perturbations 
in the longitudinal gauge. 

3.3 Modification of gravity 

Having derived the first order perturbations of the field equations along with the matter 
density perturbations, it is now possible to study how the relevant quantities in the subhorizon 
approximation deviate from the usual f(R) theories and GR. Firstly, one shows that the 
identity $ = f no longer holds: since the non-diagonal components of the energy- momentum 
tensor vanish, 5T l j = 0, the last term on the r.h.s. of the spatial-spatial component of the 
modified field equations, Eq. (3.3), must vanish for i ^ j [1]. This leads to a condition 
involving the potentials and which, for the present model, has an additional term due to the 
modification of gravity. Recalling that 5 denotes terms which are linear both in the metric 
and matter fluctuations such that 



<»' '> < h ( ^ - 2^9 ) m 2F-, S f , . (3. 10! 



one gets 



$ - = -ShxF = -<51n(Fi - 2F 2 p) , (3.11) 

which reduces to the usual GR result = ^ in the absence of a modification of gravity. 

For the formation of LSS, the relevant scales correspond to modes deep inside the horizon 
during the matter dominated era, since they reenter the horizon before matter-radiation 
equality is attained. These obey the condition k 2 /a 2 H 2 S> 1 and, noting that during the 
matter-dominated era one has R ~ H 2 , one concludes that k 2 /a 2 R S> 1. Within this regime, 
one may neglect time derivatives of the potentials and retain only terms on the energy density 
fluctuations [12, 15] (further discussion on the validity of the subhorizon approximation is 
found in Refs. [25, 36]). 

In this approximation, the fluctuation of the scalar curvature in Fourier space becomes 

k 2 

5R « -2-;r (2* - $) . (3.12) 
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For GR, the time-time component of the field equations, {k / a) 2 {2^ / k) ~ — 5p, along with 
the relation = ^ for the potentials leads to the standard generalized Poisson equation of 
the form, <3? ~ —4:irG(a/k) 2 5p, where is regarded as a generalized Newtonian potential. For 
a modified theory of gravity, one has a new coupling constant G, which differs from G due to 
the presence of factors involving the functions fi(R), as shown below. Indeed, the equation 
for the potential $ in the subhorizon approximation has the form of a Poisson equation, 



$ « -AirG^Sp , (3.13) 

where one identifies a modified gravitational coupling constant, G, which differs from G due 
to the modification of gravity. 

The expression for G is obtained from the time-time component of the modified field, 
Eq. (3.2), which, in this approximation, becomes 

- (1 + f 2 ) 5p - (F 2P ) 5R = -%[-F 1 (2$) + (2F 2 p) $ + (5F)] , (3.14) 

cr 

where 5R is given by Eq. (3.12). The fluctuation, SF, in Eq. (3.10), can be written as a 
function of the potential, <£, and the matter density fluctuation, 5p, from the perturbation 
on the scalar curvature, Eq. (3.12). Using this result in the time-time component of the field 
Eqs. (3.2), with the subhorizon approximation, Eq. (3.14), leads to an expression for the 
modified gravitational coupling constant 



G =1 + 2A(2mi _ m2 ) E 



where one identifies the quantities 



S ee I±£ = J + f * , (3.16) 



kF k(F 1 - 2F 2 p) ' 
and 

"» s 4- *fcfg ■ < 3 17 » 

m 2 = Rrr^j- ■ (3.18) 

1 + J2 

As expected, inserting fi{R) = R/k and f 2 (R) = yields m\ = m 2 = and one recovers 
the usual gravitational constant, G = G. Notice that, for f(R) theories, m 2 vanishes: this 
parameter thus emerges as a new result due to the effect of the non-minimal coupling. 

From the relation between the potentials, Eq. (3.11), the equation for the second po- 
tential has also a Poisson form, 

a 2 

* w -4irGq-^8p , (3.19) 
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and the deviation from GR is now reflected in the quantity 



1 + 2-|^ (mi + m 2 ) 

As already mentioned, the impact on weak lensing provides a straightforward test to 
these modified theories of gravity. The potential that determines the weak lensing is given by 
&WL = $ + Vl/. Using the Poisson equations for the potentials, one may rewrite Eq. (3.16) in 
the suggestive form £ = (G/G + q)/2, so that that one may write <&wl — —^G(a/k) 2 5pT 1 
[26], showing the connection between £ and a potential deviation from GR weak lensing. 

By the same token, it is also possible to construct the anisotropic parameter r\ = ($ — 
V &)/ X I / , which also indicates a deviation from GR (rj = 0), 



2-1^ (mi - 2m 2 ) 

= ggV 2J_ _ 

l + 2^(m 1 +m 2 ) 
3.4 Comparison with f(R) theories of gravity 

After obtaining the main results induced by the modification of gravity, it is now possible to 
make a comparison between the usual f\(R) theories and the non-minimal coupling with an 
additional f 2 (R) coupling function. 

First, in the equations for the modified field equations, Eqs. (3.2) and (3.3), it is clear 
that the presence of the non-minimal coupling results in a modification of the l.h.s., which 
is a characteristic of this model, and two modifications in the r.h.s. of the same equations, 
which consist in the inclusion of a new term dependent only on the coupling function and an 
extension of the term on the Fi(R) function to the linear combination, F = F\ — 2F 2 p. This 
last extension appears throughout the obtained results and, in particular, in the difference of 
the potentials in Eq. (3.11). 

The other relevant distinct property resorts to the gravitational interaction depicted 
by the X quantity, Eq. (3.16), and the parameters m\ and m 2 , given by Eqs. (3.17) and 
(3.18). These determine the form of the gravitational interaction, G, the factor, q, on the 
modified Poisson equation for the potential, in Eqs. (3.15) and (3.20), respectively, and 
the anisotropic parameter, r], in Eq. (3.21). In comparison with fi(R) theories, the coupling 
function results in a modification of mi and S, and the inclusion of a competing term m 2 . 

For the £ quantity it is clear the same extension of the F\ function to the combination 
F = F\ — 2F 2 p, and the presence of the f 2 (R) coupling function, such that, if one has a 
dominant coupling function, f 2 (R) 3> 1, then 



s - • < 3 - 22 > 

and one concludes that the presence of the coupling function overcomes the usual fi(R) 
theories of gravity in the weak lensing behaviour. 

Then, for the parameter mi, Eq. (3.17), one has also the same extension for the function 
F%, but the parameter m 2 corresponds to an emergent result in this new theory of gravity, 
and indicates the dominance of the coupling function. Hence, it is clear that the parameters 
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rrii are responsible for different regimes of dominance of the fi{R) functions, and it is S which 
in fact induces the deviation from GR in both theories. 

At last, in the context of the evolution of matter overdensities, the usual fi(R) theories of 
gravity only induce the already seen modification in the gravitational interaction, G. However, 
the inclusion of the f2{R) coupling, besides the modification of gravity as already seen, it also 
induces a new behaviour from the presence of the extra-force, which was not present in the 
former. This effect leads to the presence of the coupling potential, Eq. (3.5), and respective 
fluctuation, <5<3> c , in the equation for matter density perturbations, Eq. (3.9). In the following 
section one will see the physical implications of such terms in the cosmological dynamics of 
matter overdensities. 



4 Evolution of matter overdensities 

4.1 Equation of motion in the subhorizon approximation 

For the formation of LSS one must consider the evolution of matter density perturbations 
from the set of equations, Eqs. (3.8) and (3.9). If one introduces the gauge-invariant matter 
density perturbation, 5 m = 5p m / p m + 3Hv, the respective equation for the evolution of matter 
overdensities, in Fourier space, becomes 



k 2 r i 

S m + (2H + $> c )5 m + ^{<S> + 5® c ) = 3 B + (2H + <f> c )B , (4.1) 

where B = + Hv. This result shows that the presence of the non- minimal coupling induces 
the background term, $ c , and a sort of backreaction effective perturbative term, 5& c , due 
to the presence of the extra force. In this case, the subhorizon approximation consists in 
neglecting the terms on the time derivatives of B, so that the r.h.s. of the above expression 
may be disregarded. 

The term on S m , which is due to the expansion of the Universe in the absence of the 
coupling, has a frictional effect: its absence would lead to an exponential growth of matter 
density perturbations. In the presence of the coupling function there is a new similar term 
which shows up in the equation, <I> C = /2-R/[(l + f-z)R\ This indicates that the non-minimal 
coupling acts to reinforce (or decrease) the effect of the background expansion of the Universe. 

To obtain further insight about this term, it is relevant to consider the space component 
of the geodesic equation for a particle in the presence of the extra force, 



du a 
ds 

At zeroth order, this component reads 



+ r>'V = r • (4.2) 



™ = -(2tf + <I> c )uV , (4.3) 
as 

clearly showing the similarity between the frictional behaviours of the coupling function and 
the one from the expansion of the Universe. This behaviour accounts for the adiabatic ex- 
pansion of the Universe without an interchange of energy between matter and curvature, in 
particular in what concerns the evolution of matter overdensities. 
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To understand the origin of the second other term, 5& c , in Eq. (4.1), one must consider 
the subhorizon approximation. Since 5& c ~ 5R, and from the perturbation on the scalar 
curvature, 5R, given by Eq. (3.12), one may also identify a Poisson equation as in the case 
of the potentials and Using the equations for these potentials, Eqs. (3.13) and (3.19), 
and the modified coupling constant, G, in Eq. (3.15), one obtains 



a 2 



5<5> c = -4nG c -^8p m , (4.4) 



where the coupling constant, G c , is given by 



Gc = -f^l ( 1 + 6 ^ m2 ) E 
G = l + 3^m x 



(4.5) 



Therefore, Eq. (4.4) reveals that the extra force, through the perturbation in the coupling 
potential, acts as an effective gravitational interaction for the growth of matter density per- 
turbations. Here one observes the dependence on the quantity m^: in the absence of the 
coupling function mi = and the extra force vanishes. 

With this identification, the equation of motion for matter overdensities, Eq. (4.1), 
becomes 



b m + (2H + $ c ) 5 m - 4vr(G + G c )p m 5 m ~ . (4.6) 

From the above, one sees that the non-minimal coupling also induces a sort of backreac- 
tion gravitational effect on the evolution of matter overdensities from the action of the extra 
force. From the identification with an effective gravitational constant, G c , the perturbation 
on the coupling potential, <5$ c , due to the presence of matter density perturbations, 5p m , also 
induces an attractive interaction in the matter perturbations in a way which increases the 
growth of matter overdensities similar to the one of Newtonian gravity. 

Using the explicit expressions for the coupling constants, Eqs. (3.15) and (4.5), one has 
a total gravitational interaction, G e ff = G + G c , given by 



G 



eff 



1 + 4 



fc 2 

a 2 R 



mi — 771,2 



l + 3^m 2 



G 



l + 3^ mi 



(4.7) 



Through this expression, one may conclude that it is the quantity £ that in fact characterizes 
the modification of gravity, which is independent of the mode under consideration, and the 
expressions for mi indicate whether this is due to a modification of gravity, G, or rather 
the induced extra force, G c . Although the above splitting of the effect of the non-minimal 
coupling into a direct gravitational impact plus the action of an extra force can be regarded 
as somewhat arbitrary, it is very useful to better grasp the phenomenological implications of 
the model under scrutiny, and to physically interpret the behaviour of competing terms. 



- 10 - 



4.2 Dynamics for a power- law coupling function 

To understand the dynamics of matter perturbations in the presence of the non-minimal 
coupling, one must study the growth behaviour of matter overdensities from the respective 
equation of motion, Eq. (4.6), where the departure from the standard behaviour comes 
from the presence of the <J? C frictional term and the modification of gravity with the effective 
attractive interaction from the extra force, G c , resulting in the total gravitational coupling 
constant G e ff. 

Following previous works, one takes the specific coupling function of a power-law form, 
f2(R) = (R/R2) 71 , where R2 is a characteristic curvature scale. For this term to be relevant 
only at late times, and thus drive the accelerated expansion of the Universe, the exponent 
n has to be negative and i? 2 smaller than the current value of the scalar curvature [31]. 
Depending on the value of i? 2 , the additional terms derived from the non-minimal coupling 
may be dominant or perturbative at the time scale relevant for structure formation. These 
two options are considered in the following discussion. 

For this class of functions, the E-quantity, Eq. (3.16), becomes 



1 - 2n-jj-/ 2 

Considering a perturbative coupling function, /2<1, along with a perturbative behaviour in 
a way that the derivative terms in the trace Eq. (2.3) do not have a significant contribution, 
one may consider R ~ np and the previous expression, Eq. (4.8), can be approximated by 



I + /2 
1 - 2n/ 2 



l + (l + 2n)/ 2 ~l , 



(4.9) 



and, for a power-law coupling function with a negative exponent along with the condition, 
n < —1/2, one has E < 1. For a dominant coupling function one has instead 



1 R 

2n up 



(4.10) 



This expression, along with the restriction E > for a correct cosmological behaviour, places 
a strong condition for negative values for the exponent, n < 0, if the coupling dominates. 
For the quantities rrt\ and 7712, one has 



mi 



2n(l-n)gffi h 
1 - 2n^/ 2 



and 



'2n(l-n)^/ 2 , / 2 «1 
n-1, f 2 > 1 



(4.11) 



m 2 = n 



h 



I + /2 



n/ 2 , / 2 <1 
n, / 2 >1 



(4.12) 
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where the corresponding values, in general negative, increase from small values when the 
coupling function is perturbative, to become of order unity in the non-perturbative regime. 
For the present case, the corresponding potential can be approximated as 



i J2 Jt 

<p r = n — — Ri 

1 + f2 R 

using F2 = for a power-law function. Since the scalar curvature decreases during 

cosmic expansion, one concludes that this term acts always as a frictional term, $ c > 0, for 
a negative exponent. 

4.3 Cosmic evolution 

In general, the coupling term, responsible for the cosmic expansion, should arise only at 
late-times, such that, at the beginning of the matter dominated era, one would have a small 
deviation of GR with m, <C 1 and rriiik 2 /a?R) <C 1. As cosmic expansion evolves, not only 
does the first factor, k 2 /a 2 R, becomes approximately constant, but so do the parameters rrii. 
However, the factor k 2 /a 2 R becomes very large and mi of the order unity: therefore, one 
should consider two different regimes for the dominance of the coupling function along with 
the factors mi(k 2 /a 2 R). For a deviation of GR only due to the non-minimal coupling one 
considers fi(R) = R/n and a power-law coupling function f2(R) ~ R n - 

The equation of motion, Eq. (4.6), can be rewritten in terms of the scale factor 



nj2ji , 72 <^ 1 



n 



R 

D 1 



h » l 



(4.13) 



C + ^ (1 - LJ eff ) S' m - 127 ^ e// n m F 6 m = , (4.14) 
where a prime stands for derivative with respect the scale factor and one defines 



= ^ — r , (4.15) 

3H 2 (F 1 - 2F 2 K Pm ) ' k ; 



and an effective parameter 



4.3.1 Matter dominated era 



H' 



(4.16) 



At the beginning of the dominated matter era, for a given perturbative coupling function, one 
should have k 2 /a 2 R mi <C 1, so that the total gravitational constant can be approximated as 



G, 



eff 



G 



1 + 



a 2 R 



(mi — 47712) 



E 



l-2n{\ + n)\-f, 



(4.17) 



n/2 (with rrii <C 1). 
■~ 3 and Qm ~ 1, and 



since the parameters rrii are given by mi ~ In (1 — n) 72 and mi 

At this stage the Universe is matter dominated, with R ~ a 
the equation for matter over densities, rewritten here as a function of the number of e- folds 
N = In a, becomes 
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<h,) +l^r-^ 1 + A n(k,a)] 5 m = , (4.18) 



where 



d 2 N 2dN 2 



k 2 k 2 
A n {k, a) = ^- (mi - 4m 2 ) = -2n (1 + n) -^/ 2 ~ a 1 " 3 ™ , (4.19) 

and one has neglected 

k 2 

u eff = 2nf 2 < < 1 , (4.20) 

in the subhorizon approximation. Assuming that the density is a function of the growth 
factor, 5 ~ exp {J f g dN}, and the approximation 
solution 



df.j 



ilN 



■C /g, one has the growing mode 



-, . 6n(l + n) fc 2 r 
S~ a 5(3n-l) a lRlna j2 ; (4.21) 

which leads to the standard result, 5 ~ a, for a matter dominated Universe, in the absence 
of the coupling function, ji = 0. Therefore, the sign of the factor, 6n (1 + n) /5 (3n — 1), 
determines whether the modification of gravity, in this regime, attenuates or enhances the 
growth of matter perturbations. The result is presented in Fig. (1), where it is shown that 
the exponents in the range — 1 < n < and n > 1/3 lead to an increased growth, and for 
n < — 1 to a decrease. 

Naively, one could expect a critical behaviour in the region near n ~ 1/3. However, one 
should first notice that, for n = 1/3, the coefficient A n becomes constant and the resulting 
evolution of the fluctuation is 

5~ a l -i< l+n) ^ h , (4.22) 
which does not present such divergence issue. 
4.3.2 Curvature accelerated expansion 

At the onset of the matter dominated era one may have instead, niik 2 / 3> a 2 R, when the 
coupling function begins to dominate and the parameters rrii become of order unity. In this 
case, the total coupling constant is given by 

°eff „ °c „ 4 k<2 m l s ~ 2n k<2 1 ( 4 23) 

G G a 2 R m\ n — 1 a 2 Kp m ' 

and, with the energy density of matter at the background level, p m ~ a~ 3 , one concludes that 
the interaction linearly increases with the scale factor G e tt ~ a. 
The equation of growth, Eq. (4.14), becomes 

d 2 5 1 dS , „ , 

+ n^ + B n(k,a)5 m = , (4.24) 



dN 2 2dN 
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with 



B n (k,a) 



-6 



k 2 m\ 
a 2 R mi 



h 



6n 2 k 2 
1 — re a 2 R 



,l-3n 



(4.25) 



and, using the perturbation written as a function of the growth factor, one obtains, f 2 ~ 
B n (k,a), as long as n < 1. The growing mode solution for the matter overdensity runs with 
the scale factor as 



2 / 6n 2 fc 2 ~ 7" 

a (l-3n) lna V 1-n TJr 12 



(4.26) 



where one must have re < 1, for consistency. In this case, one sees that there is also an 
overgrowth for n < 1/3, as in the case of a perturbative coupling — otherwise, for n > 1 /3, the 
density perturbation would decrease with cosmic expansion, which corresponds to a decaying 
mode and does not describe the correct cosmological behaviour. Besides the value n = 1/3, 
there is an additional critical behaviour for the exponent re = 1, i.e. a linear coupling function, 
/2 ~ R- For this value, the parameter mi, given in Eq. (4.11), vanishes, and the gravitational 
interaction becomes 



G, 



Gr 



G G 



eff 



-12 



a 2 R 



ni2 S 



(4.27) 



6 n (1 + n) 2 n 



5{3n - 1) n - 1 




-5 - 



Figure 1. Growth behaviour induced by the presence of a non-minimal power-law coupling function 
with exponent n. The curves depict the effect of n in Eqs. (4.21) and (4.23), for the matter dominated 
era with a perturbative (dark) or dominant coupling function (dashed). 
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with G e jf > 0, as £ is negative (cf. Eq. (4.10)). However, £, the quantity which characterizes 
the weak lensing, must be positive, leading one to conclude that n = 1 does not correspond 
to a satisfactory cosmological behaviour. 

In conclusion, it is clear that the condition for a negative exponent, n < 0, is necessary 
for a correct cosmological growth of matter density perturbations in the matter dominated 
era. In fact, if the coupling term is perturbative, the growth of the density perturbations is 
less than the standard 5 ~ a by a small amount. However, when the coupling takes order 
unity values, the perturbation grows with the scale factor faster than in GR. In the opposite 
scenario, the one of a positive exponent, the dominant coupling function would decrease 
during cosmic expansion along with matter density perturbations and one would not obtain 
the observed LSS. Therefore, if one considers the modified Friedmann Eq. (2.4), such that 



n m + n c = i , (4.28) 

where one defines the curvature density parameter 

^ c = 3(fi - 2F 2Pm )lP ' (4 ' 29) 

one would expect, for this non-minimal coupling model, a transition from a matter dominated 
Universe, Q m ~ 1, into a curvature dominated one, Vt c ~ 1, thus signalling that the Universe 
enters a phase of accelerated expansion [31]. 

In what concerns the background evolution, it was shown that the observation con- 
straints on the deceleration parameter are compatible with negative values for the power-law 
exponents such as n = —4 and n = — 10 [31]. Therefore, it is expected that the cosmic evolu- 
tion of perturbations, which favor negative values for the late-time dominance of a power-law 
coupling term, would be possible within this model of background expansion. One highlights 
that there is no fine-tuning involved in the assumed non-minimal coupling driving the back- 
ground accelerated expansion of the Universe, i.e. a wide range of exponents n is compatible 
with cosmographic reconstitutions of the evolution of the Universe (as discussed in [31] and 
references within, namely [37]). 

5 Conclusions 

In this work we have obtained the linearized field equations, Eqs. (3.2) and (3.3), along with 
the equation of motion for matter overdensities in the presence of the non-minimal coupling 
between matter and geometry, Eqs. (4.1) and (4.6). 

In what concerns the cosmic evolution of density perturbations, this modified theory of 
gravity alters not only the gravitational coupling constant, but also induces a backreaction 
effective perturbation due to the presence of the extra force along with an additional frictional 
term, Eq. (4.13), which is similar to the one due to the expansion of the Universe (see Eq. 
(4.3)). This last term reveals that the non-minimal coupling modification in the background 
acts only in an indirect form and accounts for the adiabatic expansion of the Universe without 
an interchange of energy between matter and curvature, in particular, in the evolution of 
matter overdensities. 

However, as can be seen in the resulting expressions for the coupling constants G and 
G c , Eqs. (3.15) and (4.5) respectively, it is the quantity S, introduced in order to quantify 
the deviation in the weak lensing, which better captures the deviation from GR. 
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We have also studied the modification that the presence of the non-minimal coupling 
induces on the evolution of matter overdensities — in particular, for a power-law function, 
f2(R) ~ R r \ m the subhorizon approximation, k 2 <C a 2 H 2 , and for a matter dominated era 
in comparison with the standard behaviour of 5 ~ a, in Eqs. (4.22) and (4.26). 

The analysis for the regimes of a perturbative or dominant coupling function, in partic- 
ular the obtained analytical solutions for the evolution of matter overdensities, lead to the 
conclusion that the conditions for a consistent cosmological behaviour require negative values 
for the power-law exponent. On one side, for the dominant coupling function, the weak lens- 
ing condition, S > 0, along with the evolution of matter overdensities, Eq. (4.26), impose the 
conditions, n < (and n < 1/3). On the other side, the necessity of a frictional attenuation 
in the growth of density perturbations, <3? c > 0, in both regimes, Eq. (4.13), only occurs for 
negative values of the power- law exponent. 

As posited in Ref. [31], a negative exponent naturally drives a transition from a matter 
dominated Universe to a phase of accelerated expansion. Hence, this work clearly shows 
that such a mechanism is compatible with the cosmological perturbations and the ensuing 
formation of LSS. 

It should be highlighted that the non-minimal coupling does not need to be composed 
of a single power-law form, but can be considered to be a more evolved function, expressed 
as a sum of various terms, 



where R n is a characteristic curvature value. Thus, different scenarios and scales (from 
astrophysical to cosmological) should lead to a dominance of one term over the others: for 
instance, in the context of this work, the dominant term must have a characteristic curvature 
Ri of a cosmological scale. Naturally, other particular phenomena and environments can be 
described by this function, if the typical values of the curvature are such that other terms 
of the above series become relevant and induce the dynamics in each context: in particular, 
galactic dark matter can be accounted for by a power-law with exponent n = — 1 or n = —1/3, 
with the curvature scales and i?_i/3 of the order of 1/R 2 DM , with Rdm of the order of 
the radius of the dark matter halo [29]. By studying a plethora of phenomena occurring at 
distinct scales, and thus probing different dominating power-law terms of the above expansion, 
one is then able to approximate the full behaviour the non-minimal coupling function. 

A First order perturbation of the field equations 

In the following, one introduces the linearization of the modified field equations, Eqs. (2.2), 
used to obtain Eqs. (3.2) and (3.2). 

If each quantity contributes to a small deviation from the background value, the modified 
field equations, at first order, read 




(5.1) 



(l + f2)8T$ + T$6f2 



F 5R£ + R£5F- -5*}Fi 5R 

[g w V p V u - 5f}0) 5F - [W) V P V, + g^VpVv - S^U] F 



(A.l) 
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where stands for the background form of the covariant derivative and the first order 
perturbation, and the same is assumed for the d'Alembert operator □. 
For the time-time component one has 



F8R U +R U {8F) - -FrfR 



-F 



r $ + 3^ + 3H($ + 2^) + 3{H + H 2 )(2<5>) 



lc 2 



+ 3(H + H 2 )(8F) 



-6(H + 2H 2 )(2<S>) - 6H(<S> + 4tf) + 6f + 2^($ - 2f 



(A.2) 



and 



-(5 00 V Vo - D)SF - [(<5«7 00 )V Vo + (V V + □)] F = 

= -F(2$) + (2$)(F + 3HF) + ($ + 3V)F - (SF + 3H8F) + \v 2 (8F) 



3FF(2$) + (6 + 3^)F - 3H8F + -^V 2 {8F) 



(A.3) 



In the case of the spatial-spatial component, and since the background value of the 
function, F, depends only on time, the expression is given by 



F5R) 



R)8F 



l^SR 



(g tk V k Vj - Sp)5F - (Sg^VkVjF - (g lk V k Vj - %D)F 



FSRj + R)8F - ^-tiFtfR + 8UD8F + DF) - g ik V k Vj(8F) , (A.4) 



where the first term is given by 



F8R)+R)8F — -5jFi5R 



8)F 



(H + 3H 2 ) (2$) - + 6*) - * + ^V 2 ^ - -^8 ik V k Vj ($ - ¥) 

a z a z 



+ 8){H + 3H 2 ){8F) 



+ 



6(H + 2H 2 ) (2$) + 6H($ + 4*) + 6^ + -^V 2 ($ - 2tf ) 



, (A.5) 



and one has used the spatial part of the curvature tensor for the metric, Eq. (3.1), which is 
given by 



Ri 



a 2 8>, 



(H + 3H 2 ) - (2H + 6H 2 ) ($ + #)- H(<5> + 6^) - ^ + ^V 2 ^ 

a 1 



(A.6) 



These results lead to the linearized field Eqs. (3.2) and (3.3). 
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